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Recent progress in the development of a gravity-sensor-based instrument for 
determining the elevation angle of DSN antennas is described. The benefits of such 
a system include the capability to locate the Gravity Referenced Elevation Encoder 
(GREE) directly on the primary reflector (thus bypassing structural flexure and 
deformation error sources), anticipated lower maintenance costs compared to the 
present gimbai encoders , direct rep/aceabib’ty or supplementation of the present 
gimbal encoders and the utilization of off-the-shelf components to construct the 
GREE. This article includes a description of the nominal GREE design. Test results 
on a laboratory breadboard model are given. Rigid-body dynamics of the GREE are 
derived and the simulated performance in response to measured antenna vibrations 
is given. 


I. Introduction 

Large structure ground-based antennas are subject to 
pointing errors which prevent achieving the full designed 
performance of the antennas [1-3]. Major error sources 
are gravity sag, thermal gradients and hysteresis of the 
primary reflector and its supporting structure. Therefore, 
the mechanical boresight of the primary reflector is not 
directly observable from the angle encoders on the axes’ 
gimbals because of the large amount of structural sepa- 
ration. As an alternative to gimbal measurements, an in- 
strument capable of sensing the Earth’s gravity vector may 
be placed directly on the primary reflector to bypass the 
structural error sources. Knowledge of the gravity vector 
in antenna coordinates defines the elevation angle. In Sec- 
tion II, a Gravity Referenced Elevation Encoder (GREE) 
is described and the results of recently completed labora- 
tory proof-of-concept tests of noise and long term drift are 
given. The rigid body dynamics of the GREE are derived 
in Section III and the nominal controller design is given. 
Simulated performance of the GREE in a DSN disturbance 
environment is detailed in Section IV. 


II. Gravity Referenced Elevation Encoder 

The GREE is constructed of all off-the-shelf compo- 
nents. Figure 1 is a pictorial of the GREE. The major 
components are the brushless DC motor, optical encoder, 
and accelerometer. The accelerometer is mounted via an 
adapter block onto the optical encoder rotor shaft. The 
opposite end of this shaft is rotated by the motor until 
the accelerometer output is nulled. The angular rotation 
required to null the accelerometer is read on the optical en- 
coder. This angle is the elevation angle of the GREE itself. 
Knowledge of the calibration matrix between the GREE 
and the primary reflector mechanical boresight yields the 
antenna elevation angle. A block diagram of the GREE is 
shown in Fig. 2 and the laboratory unit is shown in Fig. 3. 
It is anticipated that thermal stabilization will be utilized 
in the field deployable unit to prevent temperature residual 
bias drift of the accelerometer. However, the laboratory 
model has neither heaters nor thermal insulation. 

Long term stability of the sensor bias determines the 
time between recalibrations. Test results (Fig. 4) over 27 
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days indicate a 3<r of 0.8 mdegs after detrending. A trend 
of 0.2 mdeg is apparent in the data. The suspected cause 
is cyclic microtilting of the laboratory pier. The test was 
limited to an accelerometer on hand; however a superior 
version of the accelerometer, available at a modest cost in- 
crease, is expected to reduce the long term drift by 50 per- 
cent per manufacturer's specifications. Further, the drift 
may not be linear, which means that extrapolation of the 
test results is not applicable. Although the results are en- 
couraging, longer duration tests should be completed to 
characterize the bias drift. A Hewlett-Packard spectrum 
analyzer recorded the short term noise spectrum. The rms 
noise is 0.32 mdeg at a 3.63-Hz sampling resolution and an 
accelerometer bandwidth of 250 Hz. 

III. Rigid Body Dynamics 

In this section, the GREE rotational dynamics are de- 
rived. The accelerometer is a torque-balanced pendulous 
mass sensor as shown in Fig. 5. The torque r is developed 
by internal control electronics to maintain the angular po- 
sition of the mass in accelerometer body coordinates. The 
amount of torque required is r = m g / sin 9 and the out- 
put is zero when the accelerometer sensitive axis is aligned 
with gravity. 

The GREE dynamics are derived utilizing a two degree- 
of- freedom planar model. One degree of freedom repre- 
sents the rotation axis of the optical encoder and the sec- 
ond degree of freedom is the accelerometer rotation axis. 
Figure 6 details the dynamic model where for simplicity 
the misalignments are assumed zero, i.e., the calibration 
matrix is the identity transformation. The following quan- 
tities are defined: 


Twice time differentiating the above equations yields: 


e = elevation angle 

9 1 = rotor shaft angle with respect to gravity 

Xi , yi = location of center of mass of member i, * = 

1, 2 

so, Vo = rotor axis of rotation 

CO Mi = center of mass of member t, i = 1,2 

k i = distance of accelerometer axis from rotor 
COM\ 

9] = accelerometer pendulous mass angle with re- 
spect to accelerometer case 

#2 — accelerometer pendulous mass angle with re- 
spect to gravity 

f x \) fVi — resultant coupling forces on rotational 
axis i, * = 1,2 

g = gravitational constant 

li = CO Mi offset from rotational axis i, i — 1,2 

Ii = inertia of member i about COM*, i — 1, 2 

The axis of rotation of the rotational member con- 
sisting of the optical encoder shaft, accelerometer, and 
actuator rotor (see Fig. 1) is given by the coordinates 
xo,2/o* The kinematics of the system defined in Fig. 6 
are (s, = sin#,*,^ = cos#,-): 

X\ = xo + h c x x 2 = x l + kiC\ *f I2C2 

yi = yo + 2/2 = yi + k\si + h$2 


Xi = x 0 - /, (d\ci + 

j/i = i/o + h + lid) 

(1) 

X 2 — Xq — /1 (ejci + 6*iSi^ — ki + (?lSi^ — I 2 ^2 c 2 + 6*1^2^ 

h = X 0 + /1 + 6*iCj) + ki (-0,51 + ^lC! j + / 2 + ^ 2 ^ 2 ) j 

The translational dynamics are 
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m x x x - f Xl - f X2 + m x g m x ij x = / yi - f y . 


m 2 x 2 - / X3 4- m 2 g rn 2 y 2 = fy 


(2) 


Solving for the coupling forces (/ X( , / y J in terms of the translational accelerations and substituting out these accel- 
erations via Eq. (1) yields 

f Xi = (mi + m 2 ) |io - (Q\ c i + #1*1)] — m 2 £i \o\c\ + 0iSij - m 2 /2 |^ 2 C2 4- ^2^2] - (mi 4- m 2 ) g 1 


/ Xa = m 2 Xo — m 2 /i -I- 0\8ij - m 2 ki ^4-0 2 Ci -f #iSij — m 2 / 2 p 2 c 2 4- ^ 2 ^ 2 ] “ m 2 0 

/yi = ( m i + m 2 ) t/o + (mi + m 2 ) /1 ^-0 2 si 4- 0\C\ j 4- m^ $ls x + 0iCi) -f m 2 / 2 # 2 s 2 + ^ 2 C 2 ^ 

/ ya = m 2 yo -f m 2 /i ^-0 2 s 2 + #iCi) 4* m 2 fci ^-0 2 si 4- 0iCi^ -j- mi/ 2 (-0 2 S2 + #2^2) 


(3) 


The rotational dynamics may now be written and the coupling forces substituted out to arrive at the dynamics in 
rotational coordinates. The results are 


[7i + l\ (mi + m 2 ) 4- 2l x k x m 2 4- k\m 2 ] 6\ 4- [(/1 + k x ) l 2 m 2 cos(0 2 - 6 X )] 0 2 


- [(h + ^1)^2 sin (0 2 - 0i)]02 + hs\ (m x 4- m 2 )$f + h«im 2 fii 


= +lis\ (mi + m 2 ) x Q - l x c x (mi + m 2 ) y 0 4- fci$im 2 xo - &iCim 2 t/ 0 4- n - r 2 J 


(4a) 


2 \ 


[(h + k x )l 2 m 2 cos(0 2 - 0i)]0i + [(h 4- ^i)/ 2 m 2 sin(0 2 - 0i)]8\ 


4* [h 4- m 2 / 2 ]02 4- h m 2 s 2 g = 4‘/2m 2 s 2 x 0 - / 2 m 2 c 2 j/o 4- r 2 J 


(4b) 


In the dynamics, the accelerations (xo,j/o) represent 
composite motion of the GREE due to antenna vibrations 
(£0,$)) an< * a jitter on the elevation angle (e,e). This is 
expressed as 

x 0 = Xq — /o(sin e e 4- cos e i 2 ) (5a) 

y 0 — y' Q — / 0 ( cos e e — sin 5 £ 2 ) (5b) 


parameter /o represents the assumed off-axis mounting of 
the GREE on the antenna; and the parameter l x repre- 
sents the assumed center-of-mass imbalance on the optical 
encoder rotor. These two quantities are installation and 
fabrication constraints, respectively. 

The linearized GREE model is shown in Fig. 7. The 
actuator is modeled in the frequency domain as a simple 
first-order torque generator: 


where lo is the off-axis mounting of the GREE. Table 1 
lists the parametric values utilized for the GREE. The val- 
ues are obtained from actual measurements of the GREE 
components, vendor specifications, or calculations. The 



( 6 ) 
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where k Ti k v ,r and / are the torque constant, back elec- 
tromotive force (EMF) constant, coil resistance, and coil 
inductance, respectively, and s is the complex frequency. 
The accelerometer model is a second-order critically damp- 
ed system with an appropriate scale factor ( SF ). The lin- 
earized dynamics of the rotor are derived utilizing I 2 « Ii : 

1 

e_ _ (/i+mil?) _ kr 

T 2 , ( ( s 2 +w o) 

A generalized model of bearing friction [4], Fig. 8, de- 
scribes the restoring torque as a function of rotation. This 
represents the resistance to rotation exhibited by the opti- 
cal encoder bearing action. The model is characterized by 
coulomb torque (T c ), the spatial time constant (r c ), and 
the material factor (i). All of the curves in Fig. 8 have 
equal slope (T c \r c ) at the origin of the torque-rotation 
curve and therefore an exponential approximation may be 
utilized. An analytical torque-rotation curve with hystere- 
sis is shown in Fig. 9. This model is valid for 0 > 0 and 
0 < 0. The expression for this model utilized in the sub- 
sequent simulations is 

T[fi) = sgn 0)T C + [T( 0) - sgn 0)T e ] [e~ (8) 


where T(0) is the starting torque and 0 is the angular 
change from the starting (last reversal) angle. Since a 
linearized model is required for the controller design, a 
spring constant is defined by 

AT=^A0 = ^-A0 = kA0 (9) 

dp r c 

Table 2 lists values for the GREE components. 

The linearized dynamics (Fig. 7) are utilized to com- 
pute a controller. The calculations are considerably sim- 
plified if higher order modes are neglected. Referring to 
Table 2, the motor and accelerometer dynamics are widely 
separated from the plant dynamics. Also combining the 
back EMF and friction loops into the plant dynamics yields 
the reduced-order block diagram in Fig. 10. A classical 
proportional, integral, and derivative (PID) controller is 

H = b lS + b 0 + — (10) 

s 

The rotor angle 9\ is the error signal. If, for example, 
9\{u) = 0 for all u>, then 0 = e, and the system is error 
free. The controller H must keep 9\ small. From Fig. 10, 
the closed-loop error transfer function is 


9\ kk r s 


S 3 + 

■ k T SFk r bi' 

2Cw„ + 

S 2 + 

k T SFk r bo ( 


’ k T SFk r b-i 


r 


r ”, 


V 


(ii) 


For simplicity let the closed loop response be three first 
order responses possessing break frequencies u/i,W 2 , and 
u> 3 . The controller solution is 


bi = 


r 

k t SFk r 


[uq -j- U > 2 + U>3 — 2£u> n ] 


6 0 — 


k t SFk r 


[uqu>2 + W1W3 + U>2W 3 — 


The PID controller places a zero at the origin in the 
error transfer function, Eq. (11). For the analytical mod- 
els used, a constant input will result in zero error. This is 
desirable since the nominal trajectory of a celestial track 
has a fundamental period of 1/(86,400) Hz. The error fre- 
quency response can be made flat over a particular band. 
For example, let uq = 2fl*(0.01) and u >2 = u> 3 = 27r(500) 
rad/sec. The resulting Bode diagram for the error #i(w), 
Fig. 11, is flat from 0.01 to 500 Hz with a minimum error 
attenuation of -80 dB. 

IV. Nonlinear Simulations 


b - 1 = 


k t SFk r 


[W1W2W3] 


( 12 ) 


The nonlinear equations derived in Section III are uti- 
lized to simulate the GREE response to DSN antenna 
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vibrations. The nonlinear friction model (Fig. 9) is in- 
cluded. Measured vibration data on DSS 14, as given in 
Fig. 12, are the input. Also included are the effects of 
angular jitter on the elevation angle e(t) at a 20-Hz rate 
and 0.01-rad amplitude. A fourth-order Runge-Kutta in- 
tegration at a step size of 10 msec computed the GREE 
response given in Fig. 13. The expected error is less than 
1 mdeg under these conditions. 

V. Summary 

The present state of development of an inertial instru- 
ment for DSN antennas has been described. A breadboard 
model GREE tested in the laboratory demonstrated ac- 
ceptable noise and long-term drift characteristics. Rigid 
body analytical models and simulations have been devel- 


oped for design and analysis. A nominal controller design 
maintained the simulated elevation angle error to less than 
1 mdeg when driven by measured 70-m antenna vibrations. 

The performance of a GREE under actual antenna vi- 
brations can be determined by a field test. Because the 
GREE will retrofit both mechanically and electrically onto 
existing antennas, a performance comparison study with 
the existing gimbal encoders can be designed. The recently 
completed laboratory demonstration unit can be upgraded 
for the field test by the addition of heaters, thermal insu- 
lation, mechanical mounting hardware and an electrical 
interface. Prior to the field testing, the analytical model 
needs to be upgraded to incorporate the dominant flexible 
modes and the controller upgraded to maintain the desired 
performance in the specific antenna environment. 
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Tab!* 1. Mass properties for simulation. 


Parameter 


Value 

Inertia, kg-m 2 

/, = 1.28 X 10" 2 

I 2 = 5.65 X 10~ 4 

Mass, kg 

m, = 1.02 X 10° 

m 2 = 1.20 X 10 -4 

Distance to center of mass, m 

l\ = 5.00 X 10 -3 

l 2 = 8.89 X 10 -3 

Distance to accelerometer rotation axis, m 

k\ = 1.00 X 10 -2 


Distance to elevation axis, m 

l 0 - 1.00 X 10° 


Vibration scale factor, kg-m 

k x = 5.09 X 10 -3 


Rotor transfer function, (kg-m 2 ) -1 

k r = 7.78 X 10 1 


Rotor transfer function, rad/sec 

(jjo = 1.97 x 10° 



Table 2. Component mechanical and electrical properties. 


Component 

Value 

Brushless dc motor 

Torque constant, newton-m/ampere 

k T = 1.91 X 10- 1 

Inductance, henries 

/ = 4.32 X 10~ 2 

Resistance, ohms 

r = 2.40 X 10 1 

Voltage constant, volts/rad/sec 

k v = 1.91 X 10 -1 

Half-power frequency, rad/sec 

p 3 = 5.55 X 10 2 

Accelerometer 

Scale factor, volts/rad 

SF = 2.87 x 10 2 

Break frequency, Hz 

pi ,p2 = 2.00 X 10 3 

Optical encoder friction 

Coulomb torque, newton-m 

T c - 5.30 X 10 -3 

Spatial constant, rad 

r c = 4.00 X 10“* 

Spring constant, newton-m /rad 

k - 1.32 X 10 1 











ERROR, mdeg ERROR, mdeg 




ANTENNA 

VIBRATIONS 


QUANTIZATION NOISE 




Fig. 8. Generalized force-displacement curve of solid friction 
damping. 


TORQUE (newton-m) 



Fig. 9. Torque-rotation curve for linear controller design. 


71 







ANTENNA QUANTIZATION NOISE 



Fig. 10. PID controller and reduced-order GREE model (u>jy = 32.2 rad/aec; C ~ 0.002). 



Rg. 11. Error frequency response: (a) error magnitude (smali eleva- 
tion changes) and (b) error phase (small elevation changes). 
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TIME, sec xIO 4 

Fig. 12. Measured DSS-14 (70-m) antenna vibration data for 
simulation. 
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